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ABSTRACT 

We re-examine the locus of narrow line Seyfert 1 galaxies on the Mbh^ct (black 
hole mass-bulge velocity dispersion) plane in the light of the results from large 
new optically selected samples. We find that (1) soft X-ray selected NLSls have a 
lower ratio of BH mass to (J^oiii] than broad line Seyfert 1 galaxies; this remains a 
robust statistical result contrary to recent claims otherwise; (2) optically selected 
NLSls have systematically lower Eddington luminosity ratio compared to X-ray 
selected NLSls; and (3) as a result, the locus of NLSls on the Mbr-o' plane is 
affected by selection effects. We argue that there is no single explanation for the 
origin of the Mbh^c relation; instead tracks of galaxies on the Mbh^c plane differ 
with redshift, consistent with the downsizing of AGN activity. If these results 
at face value are incorrect, then the data imply that AGNs with high Eddington 
accretion reside preferentially in relatively late type galaxies at the present epoch, 
perhaps a more interesting result and a challenge to theoretical models. 

Subject headings: galaxies: active - galaxies: nuclei-quasars: general 



1. Introduction 

How do black holes (BHs) form, how do they grow, when do they become "active" as 
in quasars and low redshift active galactic nuclei (AGNs), what is the accretion history of 
BHs, how does it relate to their active phase and how does all this relate to the forma- 
tion and evolution of galaxies? These questions have received considerable attention in the 
literature in the past five years or so, at least in part due to the discovery of the M^^-a 
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relation (Gebhardt et al. 2000a; Ferrarese & Merritt 2000; Merritt & Ferrarese 2001). The 
tight correlation between BH mass and the bulge velocity dispersion of its host galaxy im- 
plies that the ultimate fate of the two is linked. It is not obvious, however, how galaxies and 
their BHs find their way onto the Mbh^ct relation. A large number of theoretical models 
attempt to explain the observed correlation with a variety of physical processes; regulation 
of the bulge growth by the feedback from the active BH appears to be a popular one (e.g. 
King & Pounds 2003, Hopkins et al. 2005; see Mathur & Grupe 2004 for a more extensive 
list of references) . 

Do all BHs follow the same track on the Mbu-ct plane or docs it depend on BH mass, 
redshift, galaxy properties or any other parameter? It is of interest, therefore, to find loci 
of high redshift AGNs on the M^n- cr plane. Rix et al. (1999) find that z~ 2 quasars 
have higher BH mass to host galaxy stellar mass ratio compared to that at z—0; as such 
they would he above the Tremaine et al. fit (Tremaine et al. 2003) to the Mbr-ct relation. 
This implies that high redshift luminous BHs grew fast while their host galaxies were still 
in the process of assembling. At low redshift, on the other hand, galaxies presumably have 
finished growing while black holes accreting at the Eddington rate will e-fold their mass in 
a Salpeter time. Finding the locus of highly accreting BHs on the M^n-a plane is also of 
interest, therefore, to understand the origin of the Mbh-c" relation. In the local universe, a 
class of Seyfert galaxies called the narrow fine Seyfert 1 galaxies (NLSls) are known to be 
highly accreting AGNs (Pounds et al. 1995). The locus of NLSls on the Mbh-c plane is 
therefore of considerable interest and is the subject of this paper. While there may not be 
a single answer to the origin of the Mbh-c relation, e.g. it may differ with redshift, NLSls 
promise to provide at least a piece of the puzzle. 



2. NLSls and the Mbr-c relation 

Methods that work well for measuring BH masses in normal galaxies, such as gas dy- 
namics and stellar dynamics, do not work well for active galaxies; the glare of the ac- 
tive nucleus makes it difficult to use these techniques. For AGNs, reverberation mapping 
provides a powerful technique to measure BH masses and has been employed successfully 
on nearby Seyfert galaxies (Peterson 1993). Reverberation mapping, however, is time in- 
tensive, so reliable BH mass measurements have been made of only a small number of 
AGNs. Based on the reverberation mapped AGNs, Kaspi et al. determined an empirical 
relation between Mbh and the width of the H/3 emission line & optical continuum lumi- 
nosity (Kaspi et al. 2000). This simple and well calibrated relation can be easily used to 
estimate BH masses in a large number of AGNs and has been used by a number of authors 
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(McLure & Dunlop 2002; Shields et al. 2003). Measuring a poses a bigger problem, because 
the strong AGN light washes out underlying stellar absorption lines. The width of the nar- 
row [OlII] emission line is often used as a surrogate for the bulge a. It was noted by Whittle 
(1992) that the kinematics of the narrow line region of AGNs is governed by the gravitational 
potential of the host galaxy bulge, and not by the nuclear BH. Nelson & Whittle (1995a,b) 
found that FWHM([OIII]) is correlated with bulge a. The outliers in the correlation were 
sources with strong radio jets; since most AGNs are radio-quiet, without strong jets, the use 
of FWHM([OIII]) as a surrogate for bulge a appeared to be reasonable. 

Mathur et al. (2001) were the first to place NLSls on the Mbh^o" plane and found that 
they do not follow the same relation as broad line Seyfert Is (BLSls) and normal galaxies. 
They used X-ray spectral energy distribution to estimate the BH masses and [OIII] widths to 
estimate a. Using a complete sample of soft X-ray selected AGNs, Grupe & Mathur (2004, 
Paper I hereafter) confirmed the above result. Specifically, they found that for a given cr, 
NLSls have smaller BH masses compared to BLSls. Their BLSl and NLSl samples spanned 
the same range in luminosity, so were well matched. In paper I, the H/3 width was used to 
measure BH mass and again, [OIII] widths were used to estimate a. Below we consider 
limitations of, controversies about, and implications of the above result. 



2.1. Limitations 

Paper I discusses the limitations of methods to estimate Mbh and a in detail. Here 
we reiterate a few for the sake of completeness and emphasize a few more. There are 
many sources of error in using the Kaspi et al. (2000) relation to estimate BH masses: 
(1) this relation is calibrated on reverberation mapping measurements of Mbh, which itself 
is uncertain by a factor of few because of the unknown geometry of the broad emission 
hue region. (2) While there are a few NLSls in the broad line region radius - luminosity 
correlation in Peterson et al. (2000), their sample is not large enough to cover the observed 
range of Mbh in our samples. (3) Extrapolation of the Kaspi et al. relation to higher or 
lower masses would introduce an additional source of error in BH mass estimates. As a 
result, errors on individual Mbh values are large. 

The errors on individual cr values are also large for various reasons. The correlation 
between FWHM([OIII]) and a (Nelson 2000) has a large scatter which produces one source 
of error. Secondly, [OIII] lines often show some blue asymmetry, which would overestimate 
FWHM([OIII]) of the core component. We corrected for the asymmetry in the [OIII] line 
profile in Paper I; even so, the errors on a remain large. 
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Because of large errors on individual Mbh and a measurements, we have emphasized in 
Paper I that the results are statistical in nature and are robust for determining the aggregate 
properties of the samples. These, however, are the results at face-value. Confirmation of 
these results require more accurate measurements of Mbh and a. 

2.2. Controversies 

The aforementioned results, however, are controversial. While the use of FWHM(H/3) 
as a surrogate for Mbh is well accepted, the same cannot be said about FWHM([OIII]) as a 
surrogate for a. Perhaps the most important issue regarding the use of cr[oni] as a surrogate 
for a was highlighted recently by Greene & Ho (2005, GH05 hereafter), which we discuss 
below. 

Using a large sample of narrow line AGNs^ selected from the Sloan Digital Sky Survey 
(SDSS) GH05 compared the bulge velocity dispersion cr* measured using the underlying host 
galaxy spectrum and (J^ohi] using the narrow [0111] emission line. This was an important 
study as it contained a very large sample, compared to the original work of Nelson & Whittle 
(1995, 1996). Moreover, the sample selection, data reduction and analysis were performed in 
a homogeneous way. Based on this comparison, GH05 concluded that indeed, the kinematics 
of the NLR gas is dominated by the bulge gravity; the widths of the low ionization NLR 
emission hues such [SII] and [Oil] track the stellar velocity dispersion in the mean, albeit 
with substantial scatter. As sTich, they may be used as a proxy for cr*. On the other hand, 
the width of the [OIII] line is significantly broader than cr^,, so o"[oni] cannot be used as a 
proxy for o"*. However, when the blue asymmetric wing of [0111] is removed, the width of 
the core component does track o"*. This is an important conclusion because it validates the 
use of cr[oni] (after removing the blue wing) as a surrogate for cr* in previous studies and will 
help many future studies. Moreover, as mentioned above, the main result of Paper 1 (that 
NLSls have lower mass BHs than BLSls for a given a) was derived after removing the blue 
asymmetry of [OHI] lines, and is thus a statistically robust result (see below). Nonetheless, 
GH05 point out that the scatter around the cr[oni]-c* relation is large, even after removing 
the blue wing, so (7[oiii] of the core component should be used as a proxy for o"* only in a 
statistical sense, as done and emphasized in Paper 1. Boroson (2003) also pointed out the 
same using SDSS early release data, viz. the scatter in crfoni] around the Mbh-ct relation is 



^Please note that these narrow hne AGNs are not NLSls. These are AGNs in which only narrow Unes 
from the narrow lines region of AGNs are visible while the strong nuclear continuum and the broad line 
region are hidden. These are traditionally referred to as type 2 AGNs. 
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large and therefore ctjoiii] should be used as a surrogate for cr* only in a statistical sense. The 
Nelson & Whittle studies (1995, 1996) had one advantage over GH05; they could effectively 
identify the outliers in the (T[om]^(^* correlation with disturbed galaxies and/or powerful 
linear radio sources. Given the SDSS data quality and the sensitivity of the FIRST survey 
used by GH05, such identifications could not be made. One should also keep in mind that 
the GH05 study is based on type 2 AGNs; while it is reasonable to extend it to type 1 AGNs 
if orientation is the only difference between the two types, it may not be so. 

Greene & Ho (2005) then go on to find secondary drivers of the deviations of cr[oiii] from 
a*, parameterized with Aa = logo-^oiii] — logcr*. They consider host galaxy morphology, local 
environment, star formation rate, bulge velocity dispersion, radio power, AGN luminosity, 
and the ratio of bolometric to Eddington luminosity (Lboi/LEdd) as possible secondary drivers. 
They do not find any strong correlation between Act and any of these parameters (though 
note the caveat above) except Lboi/LEdd- There appears to be a mild but systematic trend 
of higher Aa in objects with higher Lboi/LEdd (formal Spearman rank correlation coefficient 
of 0.46, with a probability of chance correlation P< 0.0001) as shown in equation 3 of GH05, 
which is 

Act = (0.072 ± 0.005) log Uoi/LEdd + (0.080 ± 0.005) (1) 

for (T[oiii]=FWHM([0111])/2.35. This led GH05 to conclude that a* is overestimated in 
objects with high Lboi/LEdd; such as NLSls and high redshift quasars (Paper I, Shields et 
al. 2003). 

Given the implications of these results (§3), it is important to establish whether Ac 
is truly a function of Lbd/LEdd- GH05 have calculated Lbd using observed L[oiii] and a 
bolometric correction factor. Since their SDSS sample consists of narrow line AGNs, broad 
H/3 lines are not observed, so there is no direct handle on BH mass, and so on LEdd- GH05 use 
the Mbh-c* relation of Tremaine et al (2003) to derive Mbh and so LEdd from the observed 
values of u*. Thus a function of cr* is compared to Ao", which in itself is a function of cr* 
suggesting that the correlation between Aa and Lbd/LEdd claimed by GH05 may be a result 
of a circular argument. To investigate further whether this is indeed the case, we re- wrote 
the above correlation equation in terms of the actual observed quantities. We find that it 
translates to: 

log (J[oiii\ = A log L[07//] + S log (7* + C (2) 
where numerical values of constants A, B and C are a result of correlations between 
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Lboi and L[oiii], Mbh and cr*, and equation 1. Thus we see from equation 2 that the actual 
relations underlying equation 1 are a mild correlation between (J[oiii] and L[oiii] (with a slope 
A=0.072) and a strong correlation between crjoni] and (7*(with a slope B=0.71), rather than 
the correlation between Act and Lboi/LEdd- 

GH05 have also noted that Act correlates strongly with (T*(as can be seen from equa- 
tion 2), or with Mbh, assuming Mbr-ct relation. Since Lboi/LEdd depends upon BH mass, 
they consider whether Mbh or Lboi/LEdd is the primary driver of Aa. They conclude that 
Lboi/LEdd is the primary physical parameter because the correlation of Aa with this parame- 
ter is stronger (Spearman rank correlation coefficient = 0.46) than that with M-qu (Spearman 
rank correlation coefficient = —0.32). As shown above with equation 2, part of the correla- 
tion with Lboi/L/Edd has come about because of the correlation of o"[o///] with Lpju^ and the 
rest because of the correlation of a[oiii] with cr* (Nelson & Whittle 1996). 

Is cr[oHi] really correlated with L[oni]? A literature search showed that such a correlation 
indeed exists in the 2dF quasar sample (Corbett et al. 2003) with the probabihty of a chance 
correlation P — 0.005 (see also Whittle 1985). In fact, [OIII]A5007 is the only narrow fine 
showing a correlation with luminosity in their sample. The slope and normalization of 
the correlation are not given in Corbett et al. (2003); however, they give the slope and 
normalization of the correlation between [OIII] width and Lb, the luminosity estimated from 
the absolute photographic bj magnitude contributing to the emission line. It is interesting 
to note that the slope of their FWHM([0111])-Lb correlation is 0.12 ±0.043, consistent with 
A=0.072 in equation 2. We also looked for this correlation in the NLSl sample in Paper 1, 
using the data in Grupe et al. (2004). Indeed, the data are consistent with a mild correlation 
between crjoni] and L[ohi] with slope A=0.072. Thus it is apparent that at least part of the 
claimed result of GH05, that Act correlates with Lboi/LEdd, can be explained in terms of the 
observed correlation between (7[oni] and L[oni]- The rest is due to the correlation of a[oiir\ 
with cr* (Nelson & Whittle 1996). Note also that a strong correlation between Ljoni] and 
(7* is reported by Nelson & Whittle (1996), which may arise from the correlations between 
L[oiii] and crjoni] and between crjoni] and u*. 

We further investigated to what extent the results of Paper I would be compromised, if 
at all, if equation 1 were in fact a true correlation. There is other evidence in the literature 
suggesting that the excess [0111] width is correlated with Lboi/LEdd- Recently, Boroson (2005) 
has systematically studied the [0111] lines in a sample of 400 AGN spectra selected from the 
SDSS first data release. Unlike the GH05 sample, the spectra of the Boroson sample contain 
the broad H/3 lines, giving direct estimates of BH masses and so of LEdd- He finds that (1) 
objects with higher Eddington ratio are more likely to have large [OIII] blueshifts; and (2) 
objects with large [OIII] blueshifts have anomalously broad [OIII] emission lines. However, 
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these are not tight correlations; there are high Lboi/LEdd objects that do not show blueshifts, 
and lower Lboi/LEdd objects that do. Nonetheless, given the Boroson (2005) results, we need 
to investigate the effect of crpm]- Lboi/LEdd correlation on the result of Paper L 

As shown in Paper I, the BH masses of our two samples of BLSls and NLSls are 
significantly different. If both samples followed the Mbk-ct relation, Mbhoc cr^oiii] (Tremaine 
et al. 2003). We calculated the Mbh to C[^o"/yy] I'atio for all the objects in our sample, 
with new crpiii] values calculated using equation 1. These new values of o-piii] are used 
in Figure 1, where we plot the cumulative fraction for a K-S test of the distributions of 
logAfBH^4.02 logfjioiii] for the two populations of BLSls and NLSls of Paper I. It can be 
clearly seen that the two populations arc significantly different, with the formal K-S test 
probability of being drawn from the same population P< 0.001. A Student's t-test gives the 
probability of the two populations being similar to be P< 0.0001. We thus conclude that 
the Mbh to cr^^jj^ ratio of NLSls is statistically smaller than that of BLSls. 

Since equation 1 is defined for Act in which a[om]= FWHM([OIII])/2.35, that is what 
we used for (T[oiii] in the K-S test shown in figure 1. However, as discussed above, it is 
better to remove the contribution from the blue wing of [OIII] before measuring its width. 
In figure 2, we plot the cumulative fraction for a K-S test of logMBR— 4.02 log(7[oiii] in which 
C[oiii] is measured after removing the blue wing^. The two populations are clearly different 
with the K-S test probability of being drawn from the same population P< 0.001 (t-test 
probability P< 0.0001). GH05 do not give correlation equation for Act- Lboi/LEdd when 
c"[oiii] is measured after removing the blue wing. Nonetheless, we once again apply the 
"correction" of equation 1 to the cr[oiii] values calculated after removing the blue wing and 
perform the K-S test again. The two populations are still significantly different, with the 
probability of being drawn from the same population P= 0.003 (t-test probability P=0.0017). 

The above exercises demonstrate that the soft X-ray selected samples of NLSls and 
BLSls of Paper I are statistically significantly different in Mbh to (Jpjjjj ratio, even when 
correlation of Aa with Lboi/Lsdd, as given in equation 1, is accounted for. 



2.3. Why some NLSls lie close to the Mbh-c relation. 

As noted by many authors (Mathur et al. 2001, Ferrarese et al. 2001, Grupe & Mathur 
2004, Mathur & Grupe 2005, Barth et al. 2005), not all NLSls have relatively smaller BH 



^This is similar to Figure 4 of Paper 1, except that the statistic used is logMBH— 4.02 logcrjoiii] instead 

of logMBH- logcr[oiII] 
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mass for their bulge a. While some values of a were estimated using (T[oiii] as discussed above, 
some actual measurements of o in NLSl host galaxies also exist. Ferrarese et al. (2001) 
measured a in NLSl galaxy NGC 4051 for which BH mass is measured using reverberation 
mapping (Peterson et al. 2000). They found the source to lie close to the Mbh^o" relation. 
This, however, is no surprise because NGC 4051 lies close to the Mbh^o" relation even in 
Mathur et al. (2001) in which cxjoiii] was used as a proxy for cr*. More recently, Barth et 
al. (2005) measured a using Mg b and/or Call stellar absorption lines for a sample of 
NLSls selected from SDSS (Greene & Ho 2004). They also find that the measured a is 
not significantly different from that expected from the Mbh-ct relation. Clearly, these and 
similar such results are at odds with the expectation that NLSls have growing BHs (Mathur 
et al. 2001, Paper I) and are thus young AGNs (Mathur 2000). One may argue that all 
objects with accurate measurements of (j* lie on the M^n-a relation, so again, [0111] line 
widths must be at fault. However, the example of NGC 4051, and the fact that cT[oiii] after 
removing the blue wing does track cr*, all suggest that some other factor is likely involved in 
this apparent contradiction. 

One clue towards the reconciliation of these conflicting results comes from the work of 
Williams, Mathur & Poggc (2004). Since most large NLSl samples were soft-X-ray selected, 
they were clearly biased towards X-ray bright objects. To remedy this situation, Williams, 
Pogge & Mathur (2002) constructed a large, uniformly selected optical sample of NLSls 
from the SDSS early data release and found that only a fraction of them were detected in 
the ROSAT All Sky Survey. They performed follow up Chandra observations of ROSAT 
undetected sources and found that NLSls are a mixed bag. Not all of them are soft-X- 
ray bright or have steep X-ray spectra indicative of high values of Lboi/LEdd- Moreover, 
they found that the soft-X-ray power-law slope correlates with Lx/LEdd, with fiat spectrum 
sources having lower Lx/Lgdd (see also Grupc 2004; Lu & Yu 1999). This strongly suggests 
that not all NLSls are highly accreting sources and a large fraction of optically selected 
NLSls falls in this category. Indeed, Mathur & Grupe (2005) have shown that NLSls in 
their sample which have higher Lboi/LEdd have larger (T[oiii] than those with similar BH 
masses, but lower Lboi/LEdd, and are thus likely to be growing. 

We started to find the locus of NLSls on the M-bk-ct plane because NLSls as a class 
were thought to have large accretion rates relative to Eddington (m=M/MEdd) leading to 
large Lboi/LEdd compared to BLSls. The above results, however, show that not all NLSls 
have large rh and indicate that sample selection methods strongly influence the results. To 
investigate this further, we compared the distributions of Lboi/LEdd for three samples: the 
soft-X-ray selected sample of NLSls from Paper 1, the optically selected sample of NLSls 
from Greene & Ho (2004), and the BLSl sample from Paper I. The Lbd/LEdd values given 
in table 3 of Greene & Ho (2004) were corrected using the Mbh values from Barth et al. 
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(2004), which are better determined with higher quahty data on H/3. The Lboi/LEdd values 
for the soft-X-ray selected samples are from Grupe et al. (2004). One has to be cautious 
in comparing these samples because they do not use the same prescription to estimate Lboi- 
Greene & Ho (2004) use Lboi = 9.8AL5100 uniformly for the entire sample, while Grupe et 
al. (2004) estimate the bolo metric correction separately for each object in the sample using 
the observed spectral energy distribution. Inspection of figure 8 in Grupe (2004) reveals 
that the two bolometric corrections are consistent with each other. BH masses in all the 
three samples are estimated using H/5, so carry similar uncertainties. In figure 3 we plot the 
distributions of Lboi/LEdd for the three samples. What is seen is that the soft- X-ray selected 
NLSl sample peaks at high Lboi/LEdd; the optically selected sample peaks at lower Lboi/LEdd) 
and the BLSl sample peaks at even lower Lbd/LEdd- Indeed the mean logLboi/LEdd of the 
soft-X-ray selected NLSl sample is +0.24, that of optically selected NLSl sample is —0.45 
and that of soft X-ray selected BLSls is —0.75 (the average Lbd/LEdd of optically selected 
BLSls may be even lower). This result is consistent with Williams et al. (2004) who found 
that the mean value of soft X-ray power-law slope F decreases steadily from soft X-ray 
selected NLSls to optically selected NLSls to BLSls. 

The above results imply that soft X-ray selected NLSls are highly likely to contain 
AGNs with large Lbd/LEdd, but optically selected NLSls arc not. Since AGNs with large 
Lbd/LEdd, and so large m, are the likely ones with rapidly growing BHs, only they should 
have a distinct locus on the Mbh-ct plane, away from the Tremaine et al. relation. It should 
then be no surprise to find that NLSls with relatively smaller Lboi/Lsdd, notably optically 
selected NLSls, he close to the Tremaine et al. relation, as shown in, e.g., Barth et al. (2004). 

3. Discussion &; Conclusions 

Our result at face value, that highly accreting AGNs at low redshift have lower Mbh/c^^ 
ratios than those of AGNs with low accretion rates, has implications toward our understand- 
ing of the origin of the Mbr-c relation (the AGNs with high Lboi/LEdd are a subset of NLSls 
while BLSls have lower Lboi/LEdd)- It tells us that BHs grow rapidly in their high accretion 
mode and approach the Mbb.-c relation asymptotically. This appears to be the case at least 
at low redshift where BHs grow in well-formed bulges. 

At higher redshifts, however, the situation appears to be different. As discussed in §1, 
quasars at high redshift appear to lie above the Mbr-ct relation. This implies that the BHs 
in these quasars are likely to have already grown to their "final" mass, but their host galaxies 
have yet to grow further through interactions and mergers. Thus the tracks of high redshift 
quasars on the Mbh-c plane may be horizontal, from low to high mass bulges while that of 
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low redshift Seyfert galaxies may be vertical, from low to high BH masses. At low redshift, 
where merger rates are very low, bulges in some galaxies have grown to their "final" mass/ 
velocity dispersion, before their BHs have fully grown. This is consistent with the newly 
emerging picture of "downsizing of AGN activity" or the "anti-hierarchical BH growth" in 
which high mass BHs grow rapidly at high redshift while lower mass BHs grow at successively 
lower redshifts (Merloni 2004). Since we are probing BHs of ~ lO^M© in our NLSl sample, 
it makes sense that we catch them in their growing stage at present epoch (Mathur h Grupe 
2004). 

Thus, wc argue that there is no single answer to the origin of the Mbh^o" relation; it is 
a function of redshift. The feedback from AGN may be the primary mechanism governing 
the co-evolution BHs and galaxies at high redshift (Hopkins et al. 2005), but it may not 
work at low redshift where merger rates are exceedingly low. Of course, not all bulges are 
formed through mergers; disk/bar instability can also result in formation of pseudo-bulges 
(Kormendy & Kennicutt 2004). To our knowledge there is no theoretical model on the co- 
evolution of BHs and pseudo-bulges to explain the Afen^o" relation via feedback. Some other 
mechanism, e.g. controlled accretion rate due to capture of bulge stars by the accretion disk 
(KoUmeier & Miralda-Escude 2004), may play a dominant role at low redshift in ultimately 
placing all galaxies on the Mbh^o" relation. 

On the other hand, our result at face value may not be correct due to incorrect estimates 
of a (incorrect Mbh estimates are unlikely; Paper 1). It just may be that all galaxies, with 
dead or active BHs, with low or high accretion rates, follow the same Mbh^o" relation. 
Implications of such a result are perhaps even more interesting. As shown in Paper 1 (their 
figure 2), the distributions of BH masses of our soft X-ray selected NLSl and BLSl samples 
are significantly different, with average logMBH of NLSls being 6.9 and that of BLSls equal to 
7.9. If all these galaxies he on the Mbh-c relation, it would imply that NLSls preferentially 
reside in host galaxies with low mass/velocity dispersion bulges, or in later type galaxies. 
Since a subset of these NLSls have high m=M/MEdd, this implies that BHs with high 
Eddington ratio preferentially reside in later type galaxies at low redshifts, compared to 
AGNs with low Eddington accretion. This may be even a more challenging theoretical 
problem to address, than the origin of the Mbh^ct relation, for which there seem to have 
plenty of explanations. Some of this might simply be a gas supply issue. If the gas available 
to feed the BHs at present epoch is relatively small, the accretion rate on smaller mass BHs 
may be close to Eddington, but it would be substantially sub- Eddington for higher mass BHs. 
When we search for luminous NLSls based on their narrow H/3 widths, we are necessarily 
looking for smaller mass BHs with high accretion rates, and they are to be found in later 
type galaxies. Both the high mass and low mass BHs would be growing in mass by the same 
rate then, but the fractional growth would still be higher for the lower mass BHs. To keep 
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the high m BHs from moving away from the Mbh-ct relation, the gas supply in late type 
galaxies will have to be lower. Finding the locus of highly accreting AGNs on the Mbh-c 
plane is important either way. 

We are grateful to the anonymous referee for careful and thoughtful reading of the paper 
and useful comments. We also thank Greene & Ho for their comments. 



log M_BH ~ 4.02 log a 

Fig. 1. — Cumulative fraction of a K-S test for the distribution of logMeH— 4.02 logcT[oiii] for 
the two populations of NLSls (dashed) and BLSls (solid). The data from Grupe & Mathur 
(2004) are modified with new values of (^[oiii] calculated assuming equation 1. That the two 
populations are different remains a robust statistical result. 
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Fig. 2. — Same as Figure 1, but with (J[oiii] in which the blue wing of the [OIII] line is 
removed. The two populations are still significantly different. 
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Fig. 3. — The distributions of Lboi/LEdd for three samples: soft X-ray selected NLSls from 
Grupe & Mathur (dashed), optically selected NLSls from Barth et al. (solid), and soft X-ray 

selected BLSls from Grupc & Mathur (dotted). Given the large uncertainties in estimating 
Lboi/LEdd) the trend in the distributions of the three samples is of interest, not their exact 
values. 
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